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P  wave  is  obtained.  The  early  part  of  this  signal  Is  used  to  estimate  the  j 
parameters  of  a  vonSeggern-Bl andford  modification  of  tlie  i heor et i ca 1  Haskell 
explosion  RDP,  which  allows  apparent  discontinuities  in  source  velocity  and 
acceleration  at  the  beginning  of  the  RDP.  Synthetic  teleseismic  P  wave 
signals  are  obtained  for  this  RDP  at  a  shallow  and  deep  source  depth  with 
a  t*  of  1  sec.  Wave  form  comparisons  with  teleseismic  short  and  long  perlot 
observations  Of  events  from  Amchltka  and  Novaya  Zemlya  are  excellent.  The 
primary  theoretical  difference  between  the  comparison  events  in  the  two 
regions  appears  to  be  due  to  source  depth. 

In  order  to  determine  the  correspondence  between  generalized  ray  and 
normal  mode  theory  in  the  generation  of  Lg  phase,  we  have  begun  a  program 
of  synthesizing  seismograms  for  sources  in  simple  plane  layered  models  of 
the  lithospheric  wave  guide  using  both  generalized  rays  and  normal  modes. 

In  this  report  we  show  synthetics  for  the  fundamental  and  first  five  over¬ 
tones  of  the  Rayleigh  and  Love  normal  modes  for  a  simple  layer  over  a  half 
space  and  compare  them  with  generalized  ray  seismograms  at  various  ranges. 

A  comparison  at  128  km  for  a  far  field  trapezoid  source  history  of  5t^  = 

6t2  =  (St^  =  .2  sec  shows  good  agreement  especially  at  longer  periods. 

Because  of  the  good  agreement  at  long  periods,  we  increased  the  source 
history  duration  to  6t  =  6t2  =  dt^  =  .5  sec  and  calculated  SH  comparisons 
at  ranges  if  100  km  to'''1000  ^  using  100  km  Increments.  The  agreement  at 
all  ranges  is  extremely  good,  especially  in  the  relative  arrival  times  of 
various  body  phases  in  the  wave  train.  One  of  the  unexpected  results  was 
that  the  direct  SH  wave  received  most  of  its  contribution  from  the  Love 
fundamental  mode.  In  general,  the  higher  the  mode  the  Inter  its  contri¬ 
bution  to  the  wave  train. 

The  complete  linear  response  of  piano,  elastic  layered  solid  to  a 
shear  dislocation  is  investigated.  The  solution  is  expressed  as  a  summation 
of  generalized  rays  of  the  P,  SV  and  SU  potentials.  This  allows  the 
transient  response  to  be  obtained  upon  application  of  the  Cagnlard-de  Hoop 
technique.  Numerical  results  of  the  full  solution  containing  the  near  and 
far  field  terms  are  compared  for  a  whole  space,  half  space  and  layered 
model  with  asymptotic  solutions  to  establish  the  advantages  and  limitations 
of  approximate  methods. 

Because  of  the  presence  of  large  fault  systems  near  NTS  and  the 
Sahara  test  sites,  we  have  investigated  the  theoretical  seismic  radiation 
formed  by  an  explosion  in  the  vicinity  of  an  Infinite  fault  with  various 
slip  conditions.  The  results  show  that  it  is  possible  to  generate  S 
wave  radiation  by  detonating  an  explosive  P  wave  source  near  a  vertical 
fault  across  which  differential  slip  is  allowed  to  occur.  This  mechanism 
does  not  require  the  existence  of  regional  prestress,  but  will,  in  general, 
be  quite  inefficient  unless  the  fault  is  fairly  well  lubricated.  It  must 
be  noted  that  although  the  SH  wave  amplitude  radiation  pattern  show  a 
superficial  resemblance  to  double  couple  radiation  patterns,  the  phase  is 
quite  different  in  orientation  about  the  fault  plane. 
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I.  Summary 


The  research  completed  under  this  contract  during  the  period 
1  October  1977  -  31  March  1978  falls  within  the  two  broad  topics  of 
(1)  seismic  source  mechanisms  and  (2)  seismic  propagation  structure  and 
attenuation.  These  two  topics  are  hoth  relevant  to  the  detection  and 
discrimination  of  seismic  events. 

Using  a  deconvolution  technique,  which  makes  use  of  both  the  long  and 
short  period  observations  of  body  waves,  a  broad  band  deconvolved  explosion 
P  wave  is  obtained.  The  early  part  of  this  signal  is  used  to  estimate  the 
parameters  of  a  von  Seggern-Blandford  modification  of  the  theoretical  Haskell 
explosion  RDP,  which  allows  apparent  discontinuities  in  source  velocity  and 
acceleration  at  the  beginning  of  the  RDP.  Synthetic  teleseismic  P  wave 
signals  are  obtained  for  this  RDP  at  a  shallow  and  a  deep  source  depth  with 
a  t*  of  1  sec.  Wave  form  comparisons  with  teleseismic  short  and  long  period 
observations  of  events  from  Amchltka  and  Novaya  Zemlya  are  excellent.  The 
primary  theoretical  difference  between  the  comparison  events  in  the  two 
regions  appears  to  be  due  to  source  depth. 

In  order  to  determine  the  correspondence  between  generalized  ray  and 
normal  mode  theory  in  the  generation  of  Lg  phase,  we  have  begun  a  program 
of  synthesizing  seismograms  for  sources  in  simple  plane  layered  models  of 
the  lithospheric  wave  guide  using  both  generalized  rays  and  normal  modes. 

In  this  report  we  show  synthetics  for  the  fundamental  and  first  five  over¬ 
tones  of  the  Rayleigh  and  Love  normal  modes  for  a  simple  layer  over  a  half 
space  and  compare  them  with  generalized  ray  seismograms  at  various  ranges. 

A  comparison  at  128  kms  for  a  far  field  trapezoid  source  history  of 


6t^  =  =  5t2  =  2  shows  good  agreement  especially  at  longer  periods. 

Because  of  the  good  agreement  at  long  periods,  we  increased  the  source 
history  duration  to  dt^^  =  dt^  =  dt^  =  .5  secs  and  calculated  SH  comparisons 
at  ranges  of  100  km  to  1000  km  using  100  km  increments.  The  agreement  at 
all  ranges  is  extremely  good,  especially  in  the  relative  arrival  time.s  of 
various  body  phases  in  the  wave  train.  One  of  the  unexpected  results  was 
that  the  direct  SH  wave  received  most  of  its  contribution  from  tlie  Love 
fundamental  mode.  In  general,  the  higher  the  mode  the  later  its  contribution 
to  the  wave  train. 

The  complete  linear  response  of  plane,  elastic  layered  solid  to  a  shear 
dislocation  is  Investigated.  The  solution  is  expressed  as  a  summation  of 
generalized  rays  of  the  P,  SV  and  SH  potentials.  This  allows  the  transient 
response  to  be  obtained  upon  application  of  the  Cagniard-de  Hoop  technique. 
Numerical  results  of  the  full  solution  containing  the  near  and  far  field 
terms  are  compared  for  a  whole  space,  half  space  and  layered  model  with 
asymptotic  solutions  to  establish  the  advantages  and  limitations  of  approxi¬ 
mate  methods. 

Because  of  the  presence  of  large  fault  systems  near  NTS  and  the  Sahara 
test  sites,  we  have  investigated  the  theoretical  seismic  radiation  formed 
by  an  explosion  in  the  vicinity  of  an  infinite  fault  with  various  slip 
conditions.  The  results  show  that  it  is  possible  to  generate  S  wave  radiation 
by  detonating  an  explosive  P  wave  source  near  a  vertical  fault  across  which 
differential  slip  is  allowed  to  occur.  This  mechanism  does  not  require  the 
existence  of  regional  prestress,  but  will,  in  general,  be  quite  Inefficient 
unless  the  fault  is  fairly  well  lubricated.  It  must  be  noted  that  although 
the  SH  wave  amplitude  radiation  pattern  show  a  superficial  resemblance  to 


II.  Time  Functions  for  Nuclear  Explosions 
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Introduction 

In  recent  years  the  methods  of  computer  synthesis  of  body  waveforms 
have  provided  a  wealth  of  Information  about  earthquake  sources  and  about 
earth  structure.  The  same  methods  will  almost  certalnbe  be  equally  useful 
In  the  analysis  of  body  waves  from  bombs,  now  that  both  short  and  long 
period  data  for  large  events  has  become  available.  A  theoretical  repre¬ 
sentation  for  the  explosive  source  must  be  formulated  wlilch  Is  both  simple 
enough  to  be  used  easily  and  accurate  enough  to  predict  high  frequency 
body  waves.  A  number  of  possible  representations  have  been  proposed. 

Haskell  (1967)  devised  a  form  for  the  bomb  time  function  by  assuming  the 
boundary  conditions  at  the  source  and  fitting  some  observed  Fourier  amplitude 
spectra.  Helmberger  and  Harkrlder  (1971)  pointed  out  that  the  rise  time  was 
too  slow  In  Haskell's  model  and  importance  of  overshoot  in  the  long  period 
LRSM  data.  More  recently,  von  Seggern  and  Blandford  (1972),  used  a  simplified 
form  of  Haskell's  expression  allowing  a  faster  rise  time  and  modeled  the  Amchitka 
short  period  data  using  spectral  methods.  This  last  formulation  is  the  one 
which  will  be  discussed  here  since  it  is  the  simplest.  However,  the  free 
parameters  in  the  expression  for  the  bomb  time  function  will  be  reevaluated 
since  the  requirements  for  a  theoretical  time  function  in  a  synthetic  seismo¬ 
gram  calculation  Involving  both  short  and  long  periods  are  particularly 


% 


demanding. 
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METHODS 

The  approach  used  in  modeling  the  explosive  time  function  was  to  assume 
that  all  of  the  methods  for  modeling  earthquake  waveforms  carry  over  to  the 
explosive  source  case.  Thus,  the  synthetic  seismograms  were  represented  as 
successive  convolutions  of  an  instrument  operator,  a  Futterman  Operator 
evaluated  at  a  t*  of  1.0  and  a  free  surface  (P  +  pP)  operator.  As  usual,  the 
mantle  and  receiver  responses  were  assumed  to  be  delta  functions  for  the 
ranger  from  30°  to  90°.  The  bomb  source  time  function  was  constrained  by 
trial  and  error  modeling  of  waveform  and  amplitude  data  from  both  the 
Amchltka  and  the  Novaya  Zemlyz  test  sites.  In  order  to  eliminate  some  of 
the  difficulty  in  the  trial  and  error  process  of  finding  a  time  function, 
the  relatively  stable  Instrument  and  Q  operator  were  deconvolved  from  some 
good  quality  observed  records.  The  von  Seggern  and  Blandford  time  function 
was  used  to  fit  the  resulting  pulse.  A  full  synthetic  seismogram  calculation 
was  then  carried  out,  and  the  source  was  adjusted  until  a  reasonable  fit  to 
the  data  was  achieved. 

Analysis  and  Results 

The  deconvolution  procedure  used  to  establish  the  bomb  time  function 
convolved  with  the  free  surface  operator  is  the  one  discussed  by  Burdick 
and  Mellman  (1976) .  It  involves  using  both  a  WWSSN  short  period  record  and 
a  WWSSN  long  period  record  of  a  pulse  to  obtain  a  composite,  broad-band 
deconvolved  pulse  which  is  compatible  with  both.  Figure  1  (top)  shows  the 
deconvolved  pulse  for  a  Novaya  Zemlya  event  as  observed  at  ATL.  Several 
observed  records  from  Novaya  Zemlya  events  are  shown  next  and  at  the  bottom 
is  the  result  obtained  when  the  deconvolved  pulse  is  convolved  with  a  long 
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with  a  long  or  a  short  period  instrument.  The  good  correspondence  between 
the  observed  and  synthetic  records  demonstrates  that  the  deconvolved  source 
is  a  good  estimate  of  the  time  function-free  surface  pulse  from  Novaya  Zemly^. 

An  interesting  aspect  of  the  result  is  that  the  downswing  in  the  source 
pulse  is  slightly  larger  than  the  upswing.  This  is  difficult  to  understand 
since  the  downswing  must  be  primarily  pP  which  must  be  smaller  than  P. 

Figure  2  shows  that  this  effect  can  be  explained  if  the  source  time  function 
is  allowed  to  be  two-sided.  This  corresponds  to  allowing  B  to  be  greater 
than  zero  in  von  Seggern  and  Blandford's  expressions.  The  reduced  dis¬ 
placement  potentials  for  a  number  of  B  values  are  shown  on  the  left  of 
Figure  2.  Note  the  Increase  in  overshoot  which  corresponds  to  a  negative 
overswing  in  the  time  function  or  B  is  Increased.  On  the  right  of  Figure  2 
is  a  comparison  of  the  far  field  time  function  for  the  displacement  potentials 
with  the  deconvolved  pulse.  The  time  functions  have  been  convolved  with  a 
free  surface  operator,  so  they  should  correspond  directly  with  the  deconvolved 
pulse.  The  best  fit  occurs  for  a  B  value  of  about  2  to  5  which  corresponds  to 
a  2:1  or  3:1  overshoot  in  the  potential.  The  value  of  K  which  controls  the 
rise  time  in  the  formulation  was  set  at  5  sec 

Figure  3  compares  fully  theoretical  seismograms  computed  with  the 
von  Seggern  and  Blandford  source  with  some  additional  Novaya  Zemlya  data. 

The  synthetics  are  computed  for  the  case  of  a  deep  source  (P  -  pP  time  = 

1  sec)  and  a  shallow  source  (P  -  pP  time  -  .5  sec).  The  short  period  records 
shown  in  the  figure  were  from  smaller  bombs,  so  they  fit  the  shallow  event 
synthetics.  The  long  period  records  shown  were  from  large  bombs  and  they 
correspond  better  with  the  deep  source  synthetics.  The  Implication  is  that 
the  main  differences  in  waveshape  between  big  and  small  bombs  are  caused  by 
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a  difference  in  source  depth.  Undoubtedly  source  shape  changes  with 
event  size  as  well,  but  this  effect  is  small  compared  to  the  depth  effect. 

It  would  be  of  great  Interest  to  try  the  source  model  for  simultaneous  short 
and  long  period  records  of  a  bomb,  but  the  WWSSN  gains  are  set  so  these 
types  of  records  are  rare.  If  the  long  period  record  is  sufficiently  large, 
the  short  period  record  is  usually  off  scale. 

Figure  4  shows  a  comparison  of  data  and  synthetics  for  the  Amchitka 
blast.  Cannikin.  This  event  was  so  deep  that  the  pP  appears  as  a  visible 
feature  in  the  short  period  records.  This  feature  is  predicted  very 
accurately  by  the  synthetics.  The  long  period  records  show  a  peculiar  double 
backswlng  which  is  only  weakly  present  in  the  synthetics.  However,  the 
double  backswlng  is  caused  by  the  interaction  of  the  time  function  with  the 
near  source  velocity  structure  and  it  is  very  difficult  to  model  without 
precise  knowledge  of  this  structure.  The  important  point  is  that  again  varying 
just  the  depth  explains  most  of  the  observed  data. 

Figure  5  shows  a  comparison  of  the  observed  versus  computed  short  period 
to  long  period  amplitude  ratios  for  Cannikin  and  Novaya  Zemlya.  The  Novaya 
Zemlya  data  is  again  almost  all  from  small  events.  The  observed  ratio  is 
predicted  if  a  short  F  -  pP  time  is  used  in  the  synthetic.  If  a  longer  time 
is  used  one  predicts  the  long  period-short  period  ratio  observed  for 
Cannikin. 

Conclusions 

It  has  been  shown  that  for  appropriate  choices  of  the  parameters  in  the 
von  Seggern  and  Blandford  source  description  it  is  possible  to  fit  a  number 
of  teleseismlc  body  wave  observations.  These  include  the  long  period 
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waveform  and  the  amplitude  ratio  of  the  two.  A  much  more  significant 
conclusion  is  that  the  most  important  parameter  in  modeling  the  body  waves 
is  the  source  depth.  Since  the  depth  effect  is  so  much  more  important  than 
say  the  rise  time  or  the  overshoot  ratio,  it  may  be  the  only  reliable  para¬ 
meter  for  determining  yield. 


III.  A  Comparison  of  Genera] ized  Ray  and 


Normal  Mode  Seismograms 
David  G.  Harkrider  and  Donald  V.  Helmberger 

The  excitation  of  Lg^  waves  is  fairly  well  understood  in  the  context 
of  normal  mode  theory  for  laterally  homogeneous  media.  Because  of  known 
lateral  heterogeneities  in  the  upper  earth,  and  their  adverse  effect  on 
Lg.,  we  felt  that  it  was  important  to  determine  the  patlis  in  the  lithosphere 
wave  guide  from  the  source  to  the  receiver  for  this  energy.  Once  we  can 
predict  which  ray  paths  are  the  major  contributors  to  tlie  hg^  phases,  we 
plan  to  make  use  of  a  technique  for  generalized  rays  traveling  in  non-uniform 
wave  guides  in  order  to  evaluate  the  effect  of  varying  tlie  thickness  of  the 
lithospheric  wave  guide  between  source  and  receiver. 

In  order  to  determine  the  correspondence  between  generalized  ray  and 
normal  mode  theory  in  the  generation  of  Lg  phase  we  have  begun  a  program  of 
synthesizing  seismograms  for  sources  in  simple  plane  layered  models  of  the 
lithospheric  wave  guide  using  both  generalized  rays  and  normal  modes.  Here 
we  show  seismograms  for  a  32  km  crust  over  a  mantle  half  space  for  a  horizontal 
pure  shear  dislocation  source  at  a  depth  of  8  km  and  a  surface  receiver  at 
various  ranges.  Figures  1  and  2  are  the  synthetics  for  the  fundamental  and 
first  five  overtones  of  the  Love  and  Rayleigh  normal  modes  at  a  range  of 
128  km.  Figures  3  through  6  show  the  accumulated  sum  of  these  modes  and  the 
comparison  with  synthetics  generated  with  generalized  rays.  This  close  in 
distance  was  chosen  to  test  the  significance  of  normal  mode  theory  at  near 
distances  when  the  branch  line  contributions  are  ignored.  For  this  particular 
model  there  are  23  Rayleigh  modes  and  22  Love  modes  in  the  frequency  band, 

0  -  2  hz,  used  in  this  example. 
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Since  the  fundamental  Rayleigh  wave  dominates  the  mode  sum,  we  show  In 
Figures  7  and  8  the  accumulated  mode  sum  starting  with  the  first  higher 
mode.  Using  only  the  first  six  modes  we  demonstrate  good  agreement  between 
the  Rayleigh  wave  pulse  generated  by  both  techniques.  The  disagreement 
between  the  rest  of  the  wave  train  is  not  surprising  since  the  shear  velocity 
of  the  half  space  is  less  than  the  compresslonal  velocity  of  the  upper  layer 
and,  thus,  we  should  not  have  arrivals  corresponding  to  P  waves  in  the  mode 
sum.  The  only  contribution  of  the  P  waves  is  either  diffraction  phenomena  or 
in  the  generation  of  the  Rayleigh  pulse. 

For  Love  or  SH  waves,  the  agreement  is  good  especially  for  the  direct 
SH  wave,  l.e.  the  abrupt  spike,  and  fair  for  the  long  period  part  of  the 
later  arrivals.  Because  of  the  long  period  agreement  we  decided  to  increase 
the  source  history  trapezoid  of  6tj^  =  6t2  =  61^  =  .2  sec  to  .5  sec.  The 
resulting  comparison  between  the  generalized  ray  and  mode  sum  SH  seismograms 
for  ranges  of  100  km  to  1000  km  with  100  km  increments  are  shown  in  Figures 
9  and  10.  The  agreement  at  all  ranges  is  extremely  good  especially  in  the 
relative  arrival  times  of  various  body  phases  in  the  wave  train.  One  of 
the  unexpected  results  was  that  the  direct  SH  wave  received  most  of  Its 
contribution  from  the  Love  fundamental  mode  at  all  ranges. 

The  individual  contributions  of  each  mode  are  shown  in  Figures  11  and 
12  at  ranges  of  400  and  800  km.  The  direct  SH  wave  appears  as  a  large 
amplitude  high  frequency  arrival  preceding  the  fundamental  airy  phase.  In 
general,  the  higher  the  mode  the  later  its  contribution  to  wave  train.  In 
retrospect  this  is  not  surprising,  for  at  a  given  frequency  in  this  simple 
model  successively  higher  modes  have  higher  phase  velocities,  l.e.  smaller 
angles  of  Incidence,  and  lower  group  velocities  indicating  that  the  higher 


the  mode  the  steeper  the  angle  of  Incidence  and  the  greater  number  of 


reflections  and,  thus,  a  longer  travel  time. 

At  this  stage,  we  have  not  synthesized  generalized  rays  for  the  P-SV 
system  at  ranges  greater  than  128  km.  The  mode  sum  for  Rayleigh  modes  at 
a  range  of  1000  km  is  shown  in  Figure  13. 
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IV.  Modeling  Earthquakes  with  Generalized  Ray  Theory* 


Donald  V.  Helmberger 
and 

David  G.  Harkrider 


Abstract 

The  complete  linear  response  of  plane,  elastic  layered  solid  to  a  shear 
dislocation  is  investigated.  The  solution  is  expressed  as  a  summation  of 
generalized  rays  of  the  P,  SV,  and  SH  potentials.  This  allows  the  transient 
response  to  be  obtained  upon  application  of  the  Cagniard-de  Hoop  technique. 
Numerical  results  of  the  full  solution  containing  the  near  and  far  field 
terms  are  compared  for  a  whole  space,  half  space,  and  layered  model  with 
asymptotic  solutions  to  establish  the  advantages  and  limitations  of  approx¬ 
imate  methods. 


*Contribution  No.  2959,  Division  of  Geological  and  Planetary  Sciences, 
California  Institute  of  Technology,  Pasadena,  CA  91125 


rv-2 


Introduction 

In  the  past  decade,  significant  progress  in  our  understanding  of 
seismograms  has  been  accomplished  based  largely  on  our  ability  to  sepa¬ 
rate  source  effects  from  propagational  distortions.  Numerous  formalisms 
have  been  developed  whereby  synthetic  seismograms  can  be  computed  at 
various  positions  on  the  Earth,  for  comparisons  with  observed  seismograms. 
Iterative  techniques  can  then  be  applied  to  determine  Earth  structure 
or  source  models  or  perhaps  some  properties  of  botli.  The  most  efficient 
technique  used  to  generate  synthetics  is  essentially  controlled  by  the 
ratio  of  tiie  travel  time  of  the  disturbance  to  the  source  duration,  with 
optics  at  the  upper  limit  and  statics  at  the  lower.  However,  a  highly 
desirable  property  of  any  technique  is  tliat  the  physics  be  relatively 
apparent  so  that  the  user  can  easily  discern  how  model  parameters  influence 
the  synthetic  motion.  In  this  paper  we  will  discuss  such  a  technique, 
called  generalized  ray  theory,  where  various  useful  approximations  are 
made  based  on  the  travel  time  to  duration  ratio. 

We  assume  that  earthquakes  can  be  simulated  by  distributed  shear 
dislocations  and  that  the  Ecrth  can  be  replaced  by  a  layered  elastic 
medium;  both  assumptions  are  suspect  but  worth  consideration  as  viable 
models  at  our  present  level  of  understanding. 

First,  we  will  consider  a  shear  dislocation  in  a  whole  space  expressed 
in  cylindrical  coordinates  becasue  of  its  compatibility  with  t)ie  layiTed 
problem.  'I'he  layered  structure  complication  is  effectively  renK)ved  by 
a  generalized  ray  expansion  of  the  P,  .SV,  and  SH  displacement  potentials 
with  the  Cagniard-de  Hoop  technique  used  to  obtain  Itie  transient  response. 
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This  basic  technique  has  been  used  earlier  by  T’eker^^s  (1940)  and  his 
colleagues  and  is  sometimes  called  the  C;igniard-Tok.erl s  method  by  Ben-Henahem 
and  Vered  (1973).  Japanese  seismologists  use  a  similar  technique  and  call 
it  the  Cagniard-de  Hoop-Sato  method  (Kawasaki  et  al. .  1973).  The  basic 
idea  appears  to  have  originated  with  l^_mb  (1904)  and  has  been  modified 
for  special  purposes  by  many  researchers. 

Shear  Dislocation  Source 

Haskell  (1964)  introduced  a  shear  fault  where  a  discontinuity  in  dis¬ 
placement  across  a  fault  plane  was  allowed  resulting  in  a  double-couple 
radiation  pattern.  Following  De  Hoop's  (1958)  form  of  the  elastodynamic 
representation  Harkrider  (1976)  has  obt-ained  convenient  forms  of  displace¬ 
ments  and  displacement  potentials  for  a  number  of  different  coordinate 
systems.  The  solution  in  cylindrical  coordinates  lias  been  further  reduced 
to  a  form  suitable  for  the  application  of  Cagniard's  method  by  Helmberger 
(1974) ,  and  Langston  and  Helmberger  (1975) .  The  results  in  terms  of  the 
Laplaced  transformed  displacements  along  the  vertical,  tangential,  and 
radial  directions  are: 
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where  z,  r,  and  0  are  the  vertical,  radial,  and  polar  angle  coordinates 
respectively.  Tlie  P  wave  potential  ((Ji) ,  the  SV  wave  potential  (iO  ,  and 
the  Sll  wave  potential  (x)  are  expressed  by: 

P-wave : 
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wiicre  only  the  first  term  of  cacli  potential  is  required  in  describing  a 
pure  strike-slip  and  the  second  term  only  for  a  pure  dip-slip  orientation 
(Uarkrider ,  1976).  The  more  Important  definitions  are  as  follows: 
s  =  Laplace  transform  variable 
p  =  ray  parameter 


h  =  depth  of  source 

a  =  compressional  velocity 

6  =  shear  velocity 

p  =  density 

M  =  seismic  moment 
o 

with  tlie  orientation  constants  given  by: 

=  sin  20  cos  X  sin  ^  Y  cos  20  sin  X  sin  26 

A^CS.X.d)  =  cos  6  cos  X  cos  6  -  sin  0  sin  X  cos  26 

A3(0,X,6)  =  j  sin  X  sin  26 

A^(6,X,6)  =  cos  20  cos  X  sin  ^  ~  sin  20  sin  X  sin  26 

Aj (0 ,X, 6) =  -sin  0  cos  X  cos  6  -  cos  0  sin  X  cos  26 


whe  re 


0  =  strike  from  the  end  of  the  fault  plane 
X  =  rake  angle 
fi  =  dip  angle 


The  vertical  radiation  pattt^rns,  as  will  become  apparent  shortly. 


are  defined  by 
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The  Integrals  expressed  in  (2)  can  be  transformed  back  into  the  Lime 
domain  by  the  application  of  the  Cagniard-de  Hooi)  tcxhniciue,  see  piJJ>ejr_t 
and  Helmberger  (1972)  and  Harkrider  and  Helmberger  (1977).  For  example, 
the  field  function  defined  by 
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The  geometry  is  given  in  Figure  la  and  the  de  Hoop  contour  I'  in 
Figure  lb,  see  de  Hoop  (1960).  The  various  functions  of  p  arc  to  be 
evaluated  along  T  defined  by  choosing  those  values  of  p  which  make  i (p) 
real  and  increasing  where 

t(p)  =  pr  +  riy  Iz-hj  .  ( 


Thus , 
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Note  that  the  integrand  Is  real  Uitll 


P  = 


r  _  sln-c 
RV  "  V 


which  Is  the  ray  parameter  corresponding  to  Snell's  law. 
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’.(irt*  tV-l.  Source-receiver  geometry  and  complex  (p)  plane  with  branch 

cut  starting  at  (l/V)  and  running  out  along  the  real  (p)  axes. 

In  this  simple  case  we  have  a  closed  form  solution  for  various  values 
of  ri  since  tiie  equivalent  form  back  in  the  (u,k)  domain  lias  been  evaluated 
by  Harkrlder  (1976).  For  example, 

+  -|)j  (11) 

where  the  near-field  contribution  appears  in  terms  of  r.  However,  since 
we  need  to  evaluate  integrals  similar  to  (2)  with  complicated  complex 
integrands  later  it  should  be  noted  that  (6)  can  be  evaluated  for  various 
values  of  (t)  after  a  change  of  variable  as  proposed  by  Helmberger  (1968). 

A  relatively  fast  evaluation  of  this  type  of  integral  is  by  nonuniform 
quadrature  techniques  where  the  point  spacing  is  determined  by  the  rate 


of  change  of  the  integrand.  The  accuracy  of  such  techniques  will  be 
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discussed  later.  First  we  will  examine  some,  useful  approximations  by 

-1/2 

expanding  the  Integrand  of  (6)  in  terms  of  (t-x)  .  Note  that 
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to  first  order.  Thus,  we  can  approximate  (6)  by 
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X  n  dt  ; 


when  treating  a  high  frequency  source  with  duration,  T,  such  that 
T  <'  2pr. 

A  still  further  approximation  can  be  used  at  very  large  r's  when 
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where  t  =  R/V  and  p  =  p  and  (12)  reduces  to 
K  O 

6(t-t  ) 

C  (r,z,t)  = - -  -  (I'i) 

n  A 

called  the  first  motion  approximation.  This  approximation  is  valid  at 
teleselsmic  distances  where  the  ratio  of  travel  time  to  duration  is  0/ 
the  order  of  100  or  greater  and  has  proven  quite  useful  in  modeling  shallow 
earthquakes,  see  Langston  and  Helmberger  (1975). 

The  approximation  used  in  (12)  can  be  obtained  in  a  slightly  simpler 
way  by  returning  to  expression  (5)  and  using  the  asymptotic  form  of  the 
modified  Bessel  function,  namely 


K^(apr) 


where  =  4n^ .  Substitution  of  the  first  term  of  the  above  expression 
into  (5)  yields  (12)  after  applying  the  de  Hoop  (1960)  technique.  Note 
that  the  second  order  term  has  the  form  of  a  temporal  integration  of  the 
first  terra  or  the  form  of  a  near-field  effect.  We  will  discuss  the  high 
frequency  solution  in  the  next  section  keeping  only  the  first  term;  however, 
the  further  expansion  in  higher  order  terms  can  be  carried  out  in  similar 
fashion  and  will  be  discussed  In  the  numerical  results. 

High  Frequency  Solution  for  a  Multi-Layer  Problem 

Applying  the  method  of  generalized  reflection  and  transmission  coef¬ 
ficients  (Spencer,  1960),  it  i.s  possible  to  construct  a  representation  for 
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a  disturbance  which  has  traversed  some  layering  in  some  specified  mode  of 
propagation.  The  tangential  displacement  on  the  surface  of  a  layered  half¬ 
space  model  where  only  the  far-field  term  is  retained  becomes 


V(r,o,6,t) 


dt 
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D(t)  =  far  field  time  history 

product  of  reflection  and  transmission  coefficients  and  tlie 

summation  is  over  contributing  rays.  Numerical  evaluation  of  these 
expressions  are  discussed  in  detail  by  Helmberger  and  Malone  (1975) . 

The  high  frequency  approximations  for  the  other  components  of  motion 
are  somewhat  more  complicated  with  the  vertical  displacement  on  the  free 
surface  given  by 


W(r,o,0,t) 


4Trp  dt 
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1  =  1  \  “ 
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The  function  n^(i))  defines  the  product  of  .ill  transmission  and  reflection 
coefficients  along  the  path  from  the  source  to  tlie  receiver.  The  function 
R^^(p)  is  defined  by  Kp^(p)  or  depending  on  the  mode  of  propagation 

upon  arrival  at  the  receiver. 
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R  and  R  are  called  receiver  functions  and  are  derived  by  taking  the 
limiting  conditions  as  direct  P,  reflected  PP.  and  SP  converge  in  time  at 
a  free  surface,  see  Knopoff  et  al.  (1957)  and  Helmbe rge^  (1974). 

The  radial  displacement,  Q,  is  obtained  by  replacing  R^^  3tid  R^^  by 


Rp^  and  Rg^  defined  by 


The  reflection  and  transmission  coefficients  used  in  Ih(p)  arc  those 
defined  by  Helmbergcr  (1968).  Synthetic  responses  for  these  solutions 
will  be  discussed  after  we  examine  tlie  full  Cuy.niard  solution. 

Full  Cagn lard  Solution 

The  high  frequency  solution  discussed  in  the  last  section  has  many 
advantages  in  model  studies  due  to  its  simplicity.  However,  for  small 
values  of  (spr),  one  must  use  the  full  solution  by  applying  the  transforma¬ 
tions  used  in  deriving  expression  (6).  The  displacements  given  by  (1) 
can  be  evaluated  by  substituting  the  potentials  (2)  and  inverting  the 
various  terms  back  into  the  time  domain.  The  vertical  displacement  becomes 
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where  Wj ,  W^,  and  W3  correspond  to  a  pure  strike-slip,  dip-slip,  and  45° 
dip-slip  respectively.  The  strike-slip  response  can  be  written 
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p  for  a  whole  space. 


They  are  given  by  expression  (16)  for  a  receiver  on  the  free  surface. 
Similar  expressions  are  obtained  for  the  dip-slip  and  45°  dip-slip  cases 
with  n  =  1  and  n  =  0  respectively. 

The  tangential  displacement  is  slightly  more  complicated  because  of 
the  explicit  near-field  terms 
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for  a  receiver  on  the  free  surface.  The  dip-slip  result  is  similar  with 
n  =  1  and  where  the  factor  in  parentheses  is  reduced  by  two. 

Finally,  the  most  complicated  radial  component  is  expressed  by 


Q(r,z,e,t) 
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where 
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and 

Rpjj  =  "P  •  Rgj^  =  Hg  for  whole  space  and  given  by  expression  (18)  for 
a  receiver  on  a  free  surface.  Similar  expressions  for  the  dip-slip  result 
are  obtained  by  changing  the  n  to  1  and  the  factor  in  parentheses  reduced 
by  two , 

The  AS”  dip-slip  result  is  simply 


Q3(r,z,e,t) 
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Thus,  to  obtain  the  full  solution  requires  26  Integrations  for  each 
time  step  which  can  be  compared  with  numerical  results  obtained  from 
closed  form  whole  space  solutions  in  spherical  coordinates.  Furthermore, 
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ttif'  expressions  in  cylindrical  coordinates  are  poorly  behaved  at  small  values 
of  r  as  we  will  see  in  the  numerical  results  discussed  in  the  next  section. 

Numerical  Results 

In  this  section  we  will  present  numerical  evaluations  based  on  the 
formalism  presented  earlier  for  models  appropriate  for  shallow  earthquakes. 
However,  before  these  equations  can  be  applied  we  must  specify  the  unknown 
function,  D(t),  the  slip  history  or  the  dislocation  across  the  fault. 

Although  there  have  been  many  proposed  slip  histories  we  will  use  the  one 
suggested  by  Ohnaka  (1973)  defined  by 


D(t)  = 


l-(l+k.j,t)e 


0 


t  >  0 
t  <  0 


with  the  far  field  given  by 


D(t)  =  k|te 


t 


because  of  its  simplicity  and  its  similarity  to  the  Brune  (1970)  far-field 
source.  The  arbitrary  constant  k.j,  will  be  set  at  1  and  10  to  simulate  a 
relatively  large  earthquake,  M  >  7,  and  a  moderate  size,  M  <  6.  We  assumed 
the  seismic  moment  to  be 

M  =  4ttp  X  10^^  =  3.4  X  lO^**  ergs 
o  o 

which  is  appropriate  for  the  smaller  event  and  the  amplitudes  will  be 
expressed  in  cm.  For  earth  models,  we  chose  a  whole  space,  a  half  space, 
and  a  layer  over  a  halfspace  to  keep  the  complexity  at  a  minimum,  but 
still  test  some  commonly  held  assumptions,  see  Table  1.  That  is,  a  number 
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of  authors  have  modeled  local  earthquakes  by  reducing  the  amplitudes  of 
these  observations  by  two  and  assuming  that  the  earth  Is  a  whole  space, 
see  for  example  Trl funac  (1974). 


Thickness 

TABLE  1 

Model  Parameters 

P-Velocity 

S-Velocity 

Density 

Layer 

km 

km/s 

km/s 

g/cm^ 

1 

2. 

4.5 

2.8 

2.6 

2 

8. 

6.2 

3.5 

2.7 

As  mentioned  earlier  the  expressions  In  cylindrical  coordinates  become 
unstable  at  small  r's  since  the  near  field  P  wave  response  grows  rapidly 
and  is  mostly  cancelled  by  the  SV  wave.  This  means  the  numerical  answer 
is  the  difference  of  large  numbers,  an  unpleasant  situation.  For  this 
reason  and  to  test  the  heavy  algebra  we  computed  the  whole  space  response 
in  two  ways.  First,  by  evaluating  the  solution  in  spherical  coordinates, 
given  by  Harkrider  (1976),  and  rotating  the  vectors  into  the  directions 
r,  z,  and  0.  The  results  are  given  on  the  right  of  Figure  2,  with  the 
corresponding  solution  in  cylindrical  coordinates,  expressions  (20),  (22), 
and  (24),  given  on  the  left.  The  short  period  results  are  identical,  but 
the  larger  periods  are  somewhat  off.  Similar  synthetics  at  larger  ranges, 

A  =  16  and  32  km,  arc  Identical  for  the  short  and  long  periods  since  the 
near-field  blow  up  is  not  so  severe. 
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In  Figure  3  we  display  the  tangential  results,  V,  for  the  whole  space 
on  the  right,  half  space  In  the  middle,  and  the  asymptotic  on  the  left.  As 
might  be  expected  the  comparison  between  the  whole  space  and  half  space  is 
quite  good  using  the  factor  of  two  adjustment  although  there  is  some  distor¬ 
tion  caused  by  the  Rayleigh  pole.  On  the  other  hand,  the  as>'mptotic  solution 
is  adequate  at  the  largest  range  for  the  short  periods  since  most  Instruments 
would  not  detect  the  long  period  precursor.  The  only  significant  difference 
between  the  strike-slip  and  the  dip-slip  results  is  caused  by  the  radiation 
pattern. 

The  corresponding  comparisons  for  the  vertical  and  radial  components 
are  given  in  Figures  4  and  5.  For  these  components  the  factor  of  two, free 
surface  adjustment,  is  adequate  at  A  =  8,  but  at  larger  ranges  the  dif¬ 
fracted  P  complicates  the  situation.  On  the  other  hand,  the  comparison  between 
the  full  and  asymptotic  solutions  is  quite  good  for  the  short  periods  and 
even  the  long  periods  except  for  the  strike-slip  radial  components  at 
A  =  8,  16  km.  The  reason  for  this  is  given  by  the  strong  near-field 
dependence  for  this  particular  component  as  can  be  seen  from  expression  (24). 
The  results  for  the  45°  dip-slip  case  are  similar  to  those  in  Figure  5  with 
the  comparison  between  asymptotic  and  exact  being  slightly  better  because 
of  the  smaller  order  number  of  the  Bessel  functions. 

Actually,  it  is  somewhat  fortuitous  that  the  asymptotic  solution 
performs  as  well  as  it  does,  in  that  the  full  evaluation  of  the  far-field 
term  gives  a  poorer  approximation  to  the  full  solution.  The  reason  is  that 
the  so-called  far-field  term  contains  near-field  information  due  to  the 


unnatural  (r)  expansion.  Thus,  the  far-field  term  grows  rapidly  in  time 
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TANGENTIAL  STRIKE-SLIP 


Figure  IV-3.  Comparison  of  the  full  Cagniard  tangential  response  displayed 

on  the  left  and  middle  with  the  asymptotic  results  on  the  right. 
The  peak  amplitude  in  cm  Is  given  for  each  trace  with  the  whole 
space  results  doubled.  Model  parameters  are  specified  by  the 
half-space  in  Table  1  with  the  top  layer  removed. 


STRIKt-SLfP  VERTICAL 


STRIKt  -  SLIP  RADIAL 
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Figure  lV-4.  Vertical  and  radial  comparisons  for  the  strike-slip  orientation 
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E’igiire  IV-5.  Vertical  and  radial  comparisons  for  the  dip-slip  orientation 
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only  to  be  partially  cancelled  by  the  explicit  mar-field  tern  of  anotlu  r 
wave  type. 

The  comparison  between  the  a.symptotic  solution  and  tlie  exact  can  be 
greatly  improved  at  the  larger  ranges  by  Including  liighcr  order  terms,  for 
example;  see  Figure  6.  This  solution  does  quite  well  already  at  r  '  16  km 

with  tlie  results  for  other  orientations  being  nearly  identical  wit!i  the 
exact  after  suimning  two  terms.  On  the  other  hand,  the  short  period 
results  at  r  =  8  km  for  tin-  strike-.slip  case  are  better  than  those  pre¬ 
sented  in  Figuri;  h,  but  the  longer  period  results  are  not.  This  is  the 
normal  behavior  of  asymptotic  solutions,  but  in  this  situation  the  con- 
veigenct;  criteria  is  complicated  by  the  c.'incel  la  t  i  ons  between  tlie  various 
wave  types  as  ineiiL  toned  earlier. 

The  general  agreement  between  the  full  and  asymptotic  solution  that 
occurs  In  Figures  4  and  5  also  carries  over  into  the  layered  model  as  dis¬ 
played  In  Figure  7.  The  full  Cagniard  solution  was  generated  by  adding 
up  rays  where  the  integration  process  must  be  performed  on  each  ray 
separately.  Thus,  the  computation  time  is  much  longer  than  for  the  asymptotic 
solution  which  sums  the  rays  before  performing  the  convolution. 

Fortunately,  wi'  generally  dc>  not  need  many  rays  at  small  ranges 
because  the  internal  reflections  are  also  small  and  the  exact  formalism 
can  bi-  applied.  At  larger  ranges,  wliero  many  rays  are  required  to  model 
tiie  wave  guide  effects,  we  can  use  tlie  asymptotic  formulation. 


* 


Point  strike  -  slip  dislocation 

Tangential  Radial 


Vertical 


Figure  IV-6. 


Comparison  of  the  various  components  of  motion  for  a  strike- 
slip  orientation  at  A  =  16  km  for  a  whole  space.  The  top 
four  rows  contain  the  asymptotic  summation  after  1,  2,  4, 
and  12  terms.  The  full  solution  is  displayed  on  the  bottom. 
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FULL  CAGNIARD  ASYMPTOTIC 


Figure  IV-7. 


Comparison  of  the  full  solution  v.-ith  the  asymptotic  results 
(first  term  only)  for  a  layer  over  a  half-space  model  given 
In  Table  1. 


lV-27 


An  example  application  of  the  above  asymptotic  teclinlque  to  earth¬ 
quake  modeling  is  given  in  Figure  8.  The  observed  displacement  is  from 
a  strike-slip  earthquake  of  magnitude  5  occurring  in  the  Imperial  Valley, 
California.  This  event  was  located  in  a  region  of  known  velocity  structure 
given  at  the  top  of  the  figure.  The  only  unknown  parameters  were  the  source 
depth  and  slip  history,  D(t),  assuming  the  source  can  be  simulated  by  a 
point.  After  a  diligent  search  the  source  depth  of  6.9  km  and  a  D(t) 
specified  by  a  triangular  pulse  with  duration  of  1.5  seconds  was  found  to 
give  the  best  fit,  see  Heaton  and  Helmberger  (1977)  for  details  and  a  dis¬ 
cussion  of  the  other  components.  A  similar  study  was  conducted  on  a  larger 
earthquake,  the  Borrego  Mountain  event,  where  t’.'.o  fault  was  replaced  by  a 
distribution  of  shear  dislocations  (Heaton  and  Helmberger.  1976).  Both  of 
these  studies  involved  observations  at  considerable  distances,  30  and  60  km 
respectively,  where  the  first  term  of  the  asymptotic  solution  was  assumed 
appropriate.  However,  because  of  the  large  numbers  of  strong-motion 
instruments  presently  being  deployed  we  will  probably  obtain  close-in 
observations  of  a  large  earthquake  in  the  near  future,  and  thus  the  use¬ 
fulness  of  the  modeling  procedure  presented. 
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V.  Theoretical  and  Block  Conversion  of  Explosion 


P  Waves  into  SH  and  SV  Waves 
C.  Salvado  and  J.  B.  Minster 


For  the  purpose  of  investigating  the  possible  generation  of  S  waves 
by  purely  dllatlonal  sources  in  a  jointed  medium,  we  have  been  led  to 
consider  a  new  category  of  boundary  conditions.  The  concept  of  joint 
movement  and  block  motions  forces  us  to  abandon  the  concept  of  welded 
interfaces  with  continuous  displacements.  Since  we  are  principally  interested 
in  the  generation  of  shear  waves,  we  keep  the  normal  component  of  displacement 
continuous  across  the  Interface,  but  allow  the  tangential  component  to  be 
discontinuous.  The  design  of  a  suitable  boundary  condition  requires  that  we 
link  the  jump  in  tangential  displacement  to  the  tractions  along  the  inter¬ 
face,  which  must  themselves  be  continuous  since  no  layer  of  body  forces  is 
introduced . 


In  order  to  keep  the  problem  manageable,  and  to  derive  analytical 
solutions,  we  have  adopted  a  linear  boundary  condition  whereby  the  jump  of 
tangential  velocity  is  directly  proportional  to  the  shear  traction  across 
the  interface,  so  that 


<)>  £ 
I-*  B 


o  £  <{>  1 


where  t^  is  the  tangential  component  of  traction  uj  the  tangential  component 
of  velocity;  u,  8  are  the  rigidity  and  shear  velocity  respectively,  introduced 


here  for  homogeneity,  and  assumed  to  be  the  same  as  both  sides  of  the  inter¬ 


face  -  ^  is  termed  the  bending  parameter;  for  (t>  =  0  the  tangential  traction 


vanishes  and  the  Interface  is  perfectly  lubricated;  for  4i  =  1,  the  jump  in 


velocity  must  vanish  and  the  interface  is  perfectly  bonded  such.  Such  a 
boundary  condition  has  been  used  by  Murty  (1975,  1976). 

This  boundary  condition,  although  very  simple,  possesses  a  number 
of  desirable  properties.  It  allows  conversion  of  P  and  SV  waves  even  in 
the  case  where  the  interface  separates  identical  elastic  half  spaces; 

Murty  found  that  in  the  case  of  the  reflection  and  refraction  of  plane 
waves,  this  boundary  constitutes  an  energy  sink.  This  can  hold  physically 
only  in  the  limit  of  wavelengths  long  compared  to  the  thickness  of  the 
physical  joint,  which  explains  the  lack  of  a  characteristic  time  scale  in 
the  problem.  The  absence  of  such  a  characteristic  time,  and  thus  of  a 
characteristic  frequency  permits  us  to  treat  the  wave  propagation  problem 
using  the  Cagnlard-de  Hoop  technique. 

To  investigate  and  Illustrate  the  properties  of  such  an  interface,  we 
take  the  interface  to  be  vertical,  and  consider  a  purely  compresslonal 
(explosive)  source  at  a  short  distance  h  from  it.  We  restrict  ourselves  to 
rays  in  a  horizontal  plane  containing  the  source  (i.e.  with  a  take-off  angle 
of  90°).  Any  S  wave  polarized  in  that  plane  will  be  a  pure  SH  wave  by 
definition.  The  wave  propagation  problem  is  solved  using  a  high  frequency 
first  motion  approximation;  we  present  the  results  for  a  unit  step  source 
function,  and  for  several  values  ot  the  bonding  parameter  i)>. 

Figure  1  illustrates  the  geometry;  we  shall  ignore  for  the  moment  wave 
interactions  with  the  free  surface  so  that  the  problem  has  radial  symmetry 
about  OZ.  Figure  2  shows  the  peak  displacement  jump  as  a  function  of  radial 
distance  r,  for  various  bondings.  It  is  clear  that  values  of  ()>  much  greater 
than  .5  greatly  limit  the  relative  displacement  across  the  interface,  so 
that  we  expect  wave  conversion  to  be  most  efficient  for  low  ifi. 
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Figure  3  depicts  the  radiation  pattern  of  SH  waves  in  the  (z,  x) 
plane,  for  a  perfectly  lubricated  interface.  This  is  a  form  lobed  pattern 
which  exhibits  very  large  amplitudes  along  directions  which  correspond  to 
critical  reflections.  A  similar  phenomena  was  found  by  Burridge  al . 

(1964)  in  the  case  of  a  free  surface.  As  can  be  seen  on  Figures  4  and  5, 
the  radiation  pattern  changes  shape  with  Increasing  (Ji;  this  variation  is 
probably  too  subtle,  however,  to  be  observable  in  the  field.  Much  more 
significant,  on  the  other  hand,  is  the  large  drop  in  amplitude  as  i}i  increases. 
For  (Ji  =  .5,  the  largest  SH  amplitude  is  only  1/10  of  the  Incident  P  wave, 
this  ratio  drops  to  1/100  for  <(1  =  .9. 

The  amplitudes  plotted  on  these  figures  are  measured  perpendicular  to 
a  ray  joining  the  source  S  to  the  receiver.  Since  the  source  is  actually 
at  a  distance  h  from  the  Interface  where  wave  conversion  takes  place,  there 
exists  a  parallax  effect,  and  thus  a  component  of  "radial"  motion  associated 
with  the  converted  S  wave.  This  component  is  plotted  as  a  dotted  line  on 
the  figures. 

It  appears,  therefore,  that  it  is  possible  to  generate  SH  radiation  in 
the  earth  by  detonating  an  explosive  source  near  a  vertical  joint  across 
which  differential  slip  is  allowed  to  occur.  This  mechanism  does  not  require 
the  existence  of  regional  prestress,  but  will,  in  general,  be  quite  inefficient 
unless  the  joint  is  fairly  well  lubricated.  It  must  be  noted  that  although 
the  SH  amplitude  radiation  pattern  show  a  certain  resemblance  with  double 
couple  radiation  patterns,  the  phase  is  quite  different.  In  these  figures 
the  phases  are  antisymmetric  about  the  interface,  and  are  symmetric  about 
the  X  axis,  perpendicular  to  the  interface.  Further  inferences  about  S  wave 
generation  using  such  a  mechanism  must  await  further  calculations;  In 
particular,  crucial  problems  are  the  generation  of  Love  waves,  and  the 
behavior  of  such  Interfaces  limited  in  space. 
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